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ABSTRACT 

I present results from the modeling of stellar bars in nearly 300 barred galaxies in 
the local universe through parametric multi-component multi-band image fitting. The 
\ surface brightness radial profile of bars is described using a Sersic function, and param- 

eters such as bar effective radius, ellipticity, boxiness, length and mass, and bar-to-total 
luminosity and mass ratios, are determined, which is unprecedented for a sample of 
this size. The properties of bars in galaxies with classical bulges and pseudo-bulges 
O ' are compared. For a fixed bar-to-total mass ratio, pseudo-bulges are on average signif- 

icantly less massive than classical bulges, indicating that, if pseudo-bulges are formed 
t/^ . through bars, further processes are necessary to build a classical bulge. I find a corre- 

lation between bar ellipticity and boxiness, and define a new parameter as the product 
of these two quantities. I also find correlations between this product and normalised 
^SJ ■ bar size, between the sizes of bars and bulges, and between normalised bar size and 

^> ' bulge-to-total ratio. Bars with different ellipticities follow parallel lines in the latter 

0^ , two correlations. These correlations can arise if, starting off with different normalised 

sizes and ellipticities, bars grow longer and stronger with dynamical age, as a result 
I of angular momentum exchange from the inner to the outer parts of galaxies, consis- 

tent with previous theoretical predictions. A plausible consequence is that bar pattern 
CO ' speeds should become lower with bar dynamical age, and towards galaxies with more 

\ prominent bulges. 
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1 INTRODUCTION evidence that barred galaxies show flatter chemical abun- 

dance (O/H) radial gradients (Zaritsky et al. 1994, further. 
Many recent studies, horn observational and theoretical nj- o r> i r,r, ^ c i j_ ^ ii i ^-i n 

Martin & Koy 1994 nnd that the stronger the bar the flatter 
viewpoints, have established that stellar bars m disc galaxies i- i\ i i ■ i ^ i ^ c ^ 

° the gradient) and higher central concentrations oi molecular 

can play an important role m galaxy evolution see bellwood „ , ^ ^ , ir,r.r,\ rm ■ ^ r ^ 

n . T, H OT, . ^ , . gas (CO - Sakamoto et al. 1999). This movement of gas to 

& Wilkinson 1993: Kormendy & Kennicutt 2004; Gadotti f, , . , , • -u j , • 

„„„„ ^ . , . , , . , . , , the centre might m principle help build a young and kme- 

2009a, tor reviews), theoretical work indicates that the re- i- n i i j_ n i i ^ • j- i-i i 

' matically cold stellar bulge component, i.e. a disc-hke bulge 

distribution oi angular momentum, induced by the bar, m / a^i i r,r,r,r\ t j i i j_- ^ j_i ^ 

(see Athanassoula 2005). Indeed, observations suggest that 
the galaxy interstellar medium, as well as m the stellar and i. ,.i i , . , , i r , ■ i- i i , 

disc-hke bulges exist and have formation processes linked to 
dark matter components, has a number of important conse- , • i i- • ^ i -i-^- i i i j_ 

, n HI I , dynamical disc instabilities, such as bars, as opposed to the 

quences (e.g. Athanassoufa & Misinotis 2002; Athanassoufa , , , , . , • i, i , i • i i i / „ , , 

^ old and kmematically hot cfassicaf bufges (e.g. Caroffo et af. 

2003; Martinez- Vafpuesta et al. 2006; Berentzen et al. 2006). o ^ \ ■ r,r,r,i n • ^ i r,^no ti- i 

^ , , , , , , , , 1997; Gadotti & dos Anjos 2001; Erwm et al. 2003; Fisher 

Gas lying beyond the bar ends is driven outwards, whereas nnn/- t~v c i nnn-? i c onno n j ii- 

, , , , , , , , 2006; Drory & Fisher 2007; Fisher & Drory 2008; Gadotti 

gas lying withm the bar ends is driven to the centraf re- nr^^n, i r .i • \ 

, ° , ^ , ^ . i„o,- A , 2009b, and references therein), 

gions (e.g. bchwarz 1981; Combes & Germ 1985; Athanas- 
soula 1992; Friedli & Benz 1993; Friedli et al. 1994; Finer Theory also suggests how bars evolve with time, 
et al. 1995). This secular evolution scenario has been par- Broadly speaking, bars slow down their pattern rotation 
tially conflrmed, at least qualitatively, with observational speed, and get longer and thinner (i.e. more eccentric and 

stronger) during the course of their evolution, capturing 
stars from the disc (see Athanassoula 2003). Observations 
* E-mail: dgadotti@eso.org suggest that the strong bar in NGC 4608 has increased in 
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mass by a factor of « 1.7, through the capture of ~ 13% of 
the disc stars (Gadotti 2008). In addition, more evolved bars 
also show more rectangular-like face-on isophotal shapes, 
i.e. they are more boxy (as in the iV-body simulations of 
Athanassoula & Misiriotis 2002). In detail, however, sim- 
ulated bars can become abruptly shorter and thicker a few 
Giga-years after their formation, due to the onset of dynam- 
ical vertical instabihties that originate box/peanut bulges. 
(These seem to be simply the inner parts of bars that buckle 
off the disc plane and can be seen in iiiclincxl systems.) 
Then, about 1 Gyr later, they recover the original evolu- 
tionary trends. EXirthermore, a substantial gas component 
in the disc can also complicate the picture of the evolution 
of bar properties, halting how these properties change, and 
in some cases even reversing the trends (e.g. Bournaud & 
Combes 2002; Bournaud et al. 2005; Debattista et al. 2006; 
Berentzen et al. 2007). To date, there is no study aimed di- 
rectly at providing observational evidence on how bar prop- 
erties change with time. This is partially due to the difficulty 
of estimating bar dynamical ages. Although some first steps 
have been done in this direction (see Gadotti & de Souza 
2005, 2006; Perez et al. 2009), the results are as yet incon- 
clusive, and the methods developed require large amounts 
of telescope time. 

Moreover, bars arc found very often in disc galaxies, and 
the fraction of disc galaxies hosting bars seems to increase 
with time, i.e. the fraction is lower at redshift 2 ~ 1, as com- 
pared to z ^ (see Sheth et al. 2008, and references therein, 
but see also Barazza et al. 2008 and Jogcc ct al. 2004). In ad- 
dition, weaker, but also global non-axisymmetric structures, 
such as oval distortions in the disc, can as well efflciently pro- 
duce such redistribution of angular momentum. Therefore, 
the consequences of the presence of such structures should 
be conspicuous, and studies on the properties of bars and 
their host galaxies at « 0, as well as higher redshifts, can 
give direct clues on galaxy evolution. 

The first and easier step in observational studies of 
barred galaxies is to identify bars, and there are a num- 
ber of studies on the fraction of disc galaxies with bars (e.g. 
Eskridge et al. 2000; Barazza et al. 2008; Marinova et al. 
2009). Other studies have used ellipse fits to the images of 
barred galaxies to obtain bar properties such as ellipticity 
and length (e.g. Marinova & Jogee 2007; Barazza et al. 2008, 
2009), although it has been shown that ellipse fits can lead 
to an underestimation of the bar ellipticity (or an overes- 
timation of the bar axial ratio, Gadotti 2008). In Martin 
(1995), bar axial ratios and lengths were visually assessed, 
and a relation was found between bar length and the nor- 
malised diameter of the bulge. Ten years later, Erwin (2005) 
measured bar lengths using ellipse fits and found that bar 
size scales with disc size, and confirmed previous results that 
bars in early-type disc galaxies are clearly larger than those 
in late- type galaxies (see also Agucrri ct al. 2009). He went 
further and argued that this observational evidence can be 
qualitatively consistent with the simulations that show bars 
getting larger with time, if indeed secular evolution produces 
early-type disc galaxies from late-type ones (see also Priedli 
& Benz 1995; Martin & Fricdli 1997; Gadotti & de Souza 
2005, 2006). Applying ellipse fits to 2MASS images of 151 
spiral galaxies, Menendez-Delmestre et al. (2007) obtained 
measurements of bar length and axial ratio, and found a 



weak trend of higher ellipticities for larger bars, which is 
also consistent with the results from simulations. 

More recently, Durbala et al. (2008, see also de Jong 
1996; Laurikainen et al. 2005, 2007; Gadotti 2008) used 2D 
bulge/bar/disc decompositions of 97 Sb, Sbc and Sc galax- 
ies, employing i-band images from the Sloan Digital Sky 
Survey (SDSS) to obtain parameters such as bar length and 
bar-to-total luminosity ratio. Weinzirl et al. (2009) did a 
similar work using _ff-band images of 143 spirals and ex- 
plored also the stellar mass content in bars. About 60% of 
the galaxies in both samples are barred. In this paper, 1 
explore the results of 2D bulge/bar/disc decompositions of 
291 barred galaxies, using SDSS images in the g and i bands, 
from Gadotti (2009b, hereafter Paper 1). This allows me to 
study bar properties in a level of detail which is unprece- 
dented for a sample of this size. A thorough characterization 
of bars in massive galaxies in the local universe is thus put 
forth in Sect. 3.1, after a description of the data at hand 
in the next section. In Sect. 3.2, the properties of bars in 
galaxies with classical and pseudo-bulges are compared. In 
Sect. 3.3, I explore correlations between bar properties in 
order to test the predictions from simulations on the secular 
growth of bars. These results are discussed in Sect. 4, while 
Sect. 5 summarises the paper. 



2 DATA 

The reader is referred to Paper 1 for a detailed account of 
the sample selection and image decomposition. Here I sum- 
marise the most relevant aspects of these procedures to the 
present study. 

In Paper I, I have performed careful and detailed image 
fitting of all galaxies in a sample of 946 systems, firom bulge- 
less to elliptical galaxies. The sample was designed to be con- 
comitantly suitable for structural analysis based on image 
decomposition and a fair representation of the galaxy pop- 
ulation in the local universe. It was drawn from all objects 
spectroscopically classified as galaxies in the SDSS Data Re- 
lease Two (DR2) at redshifts 0.02 ^ 2 < 0.07, and with stel- 
lar masses larger than 10^" Mq. This parent sample is thus 
a volume-limited sample of massive galaxies, i.e. a sample 
which includes all galaxies more massive than 10^° Mq in the 
volume defined by the redshift cuts and the DR2 footprint. 
In order to produce reliable decompositions, and avoid dust 
and projection effects, I have applied another important se- 
lection criterion to produce the final sample: it contains only 
galaxies close to face-on, i.e. with an axial ratio b/a > 0.9, 
where a and h arc, respectively, the semi-major and semi- 
minor axes of the galaxy at the 25 g-band mag arcsec"^ 
isophote. This criterion also eases the identification of bars, 
which are difficult to see in very inclined projections. I have 
found that the final sample is representative of the local 
population of massive galaxies. This was done by compar- 
ing the distributions of several main galajcy properties, such 
as absolute magnitude, D„(4000)^ and concentration, in the 

^ The Dn (4000) index is based on the 4000A discontinuity seen 
in optical spectra of galaxies, and it is not sensitive to dust at- 
tenuation effects. Young stellar populations have low values of 
D„(4000), as compared to old stellar populations. 
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volume-limited and final samples, and verifying that these 
distributions are similar. 

Two-dimensional fits were performed using the BUDDA 
code (de Souza et al. 2004; Gadotti 2008) and SDSS images 
in the g, r and i bands, including up to three components 
in the models, namely bulge, disc and bar.^ The presence 
of these components was assessed by individual inspection 
of images, surface brightness radial profiles and isophotal 
maps. Given the typical redshift of the galaxies in the sam- 
ple, z ^ 0.05, and the typical FWHM of the PSF in SDSS 
images, FWHM r; 1.5 arcsec, the typical physical spatial 
resolution in these images is thus 1.5 kpc. This means that 
I likely missed most bars with semi-major axis shorter than 
Lbar ~ 2—3 kpc, typically seen in very late-type spirals (later 
than Sc - Elmegreen & Elmegrecn 1985). The whole of these 
faint bars is typically within 2 — 4 seeing elements, and thus 
they do not imprint clear signatures in either the isophotal 
contours or the intensity profile. The results presented here 
thus concern the typical, bonafide bars seen in early-type 
spirals and lenticulars. Note that, in contrast, bulges with ef- 
fective radius re of the order of one PSF HWHM can still be 
identified in the intensity profile, since they usually contain 
a much larger fraction of the galaxy light than these short 
and faint bars. The final sample includes 291 barred galaxies, 
which are the subject of the present paper. Classical bulges 
and pseudo-bulges are separated using the Kormendy (1977) 
relation, where pseudo-bulges can be identified as outliers in 
an objective fashion. It is worth noting that pseudo-bulges 
here refer to disc-like bulges, not box/peanut bulges, since 
the sample contains only face-on galaxies, and bars are in- 
cluded in the fitted models (see Paper I). 

Bulge, disc and bar are described as concentric ellipses, 
which can have different position angles and cUipticities. 
They also follow a surface brightness radial profile: in the 
case of the disc, this is an exponential (Freeman 1970) pro- 
file. Bulge and bar follow a Scrsic (1968) profile, in which the 
Scrsic index is a free parameter that controls the shape of 
the profile. A large Sersic index corresponds to a highly cen- 
trally concentrated profile with important wings, whereas a 
small Sersic index corresponds to a flatter profile with a fast 
declining outer part. It should be noted that, although both 
bulge and bar have surface brightness profiles following the 
same parametric function, there is little space for degener- 
ate solutions, since bars are generally more extended and 
eccentric than bulges, and have lower Sersic indices. 

Since I have done multi-band decompositions, I was able 
to estimate the g — i integrated colour of each component 
separately. Using the relation between g — i and the stel- 
lar mass-to-light ratio in the i-band from Kauffmann et al. 
(2007), I have all parameters necessary to accurately calcu- 
late the stellar masses of all components, including the bars, 
of the galaxies in the sample. 

The structural parameters obtained in Paper I used in 
this work are thus: 

• disc scale length /i; 

• bulge effective radius re (i.e. the radius that contains 
half of the light coming from the bulge); 

• bulge-to-total luminosity ratio B/T; 

^ Results are available at http://www.mpa- 
garching.mpg.de/~dimitri/buddaonsdss/buddaonsdss.html. 



• bar effective radius re,bar (i.e. the radius that contains 
half of the light coming from the bar); 

• bar-to-total luminosity ratio Bar /T (lum.); 

• bar-to-total mass ratio Bar/T (mass); 

• bar Sersic index near; 

• bar ellipticity e, defined as 1 — b/a, where a and b are, 
respectively, the bar semi-major and semi-minor axes; 

• bar boxiness c, defined through the equation of the gen- 
eralised ellipse used to fit the bar (see Athanassoula et al. 
1990): 



where, again, a and b are, respectively, the bar semi-major 
and semi-minor axes, and x and y are the pixel coordinates 

of the ellipse points; 

• bar length (semi-major axis) Lbar. 

It is worth noting that the bar length as determined through 
the image decompositions is usually larger than the radius of 
the peak in ellipticity inside the bar, a parameter commonly 
used to define bar size, estimated through ellipse fitting. 
However, Lbar agrees well with two other parameters used 
to define bar size, namely the first minimum in ellipticity 
outside the ellipticity peak in the bar, and the point where 
the position angle of the fitted ellipses differ by more than 
10° from the position angle of the bar (see Gadotti 2008, 
fig. 4), which are known to provide more accurate estimates 
(see discussion in Erwin 2005). 

Two other parameters used in this paper were obtained 
from the SDSS database: 

• r24- the radius of the galaxy 24 r-band mag arcsec"^ 

isophote; 

• R90: the radius containing 90 per cent of the total 

galaxy light. 

All these parameters refer to measurements using the i- 
band image, unless otherwise noted. In the remaining of the 
paper, I will explore these physical properties firstly to pro- 
vide a comprehensive description of the structure of stellar 
bars in massive galaxies in the local universe, and, secondly, 
to verify if one can test the theoretical predictions described 
in the previous section about how such bars evolve in time. 



3 RESULTS 

3.1 The structural properties of bars 

As discussed at the Introduction, large-scale stellar bars play 
a fundamental role in galaxy evolution. It is thus clearly 
very important to put forth a description of their structural 
properties as detailed as possible. With the data described 
above, I do this in Fig. 1. This gives us a thorough portrayal 
of bars in the local universe, which is useful in many ways. 
For instance, if one wants to put bars ad hoc in a theoreti- 
cal framework to describe barred galaxies today, it is most 
likely that one would want such bars to be consistent with 
those studied here. Conversely, theoretical studies aiming at 
the evolution of bars in time should be able to explain the 
existence of bars today with the properties depicted in Fig. 
1. In addition, a comparison between the properties of bars 
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Figure 1. Distributions of several bar properties. Prom top to bottom and left to right: Sersic index, ellipticity, boxiness, length (semi- 
major axis), length normalised by disc scale-length, length normalised by the radius of the 24 r-band mag arcsec"^ isophote, length 
normalised by the radius containing 90 per cent of the total galaxy light, bar-to-total luminosity ratio, and bar-to-total mass ratio. 
Marked at each panel are the median and standard deviation values of the corresponding distribution, as well as the mean la error of a 
single measurement, when available. Bin sizes are ss 1 — 2a. 



at different redshifts gives us directly a way to see how bars 
evolve in time. 

Some features in Fig. 1 are worth mentioning. First, the 
typical value of the bar Sorsic index (near ~ 0.7) depicts a 
profile which is about half way between a Gaussian (Sersic 
index equals 0.5) and an exponential (Sersic index equals 1). 
The distribution of bar ellipticities peaks at « 0.6, which is 
about 20% higher than the peaks found in studies based on 
ellipse fits (e.g. Menendez-Delmestre et al. 2007; Marinova & 
Jogee 2007; Barazza et al. 2008; Marinova et al. 2009). This 
difference is due to the fact that the isophotes enclosing bars 
are slightly rounder than the bar itself, due to the contribu- 
tion of the axisymmetric light distribution from bulge and 
disc, and thus the ellipse fits to those isophotes will also be 
rounder than the bax. In 2D fits, however, the light from 
bulge and disc is taken into account and thus a better esti- 
mate of the bar ellipticity can be obtained in this way. This 
effect was clearly demonstrated in Gadotti (2008). In fact, 
it was shown that ellipse fits indicate bar ellipticities which 
are, on average, 20% lower than those measured using 2D 
fits, employing the same galaxy images. Furthermore, Mari- 
nova et al. (2009) found that their ellipticity estimates are 
higher on average in bulgeless galaxies. Also, most bars are 
not well fitted with a pure ellipse, i.e. a generalised ellipse 
with c = 2, since bar boxiness peaks at c = 3. This is im- 



portant to keep in mind when one does 2D fitting of images 
of barred galaxies. 

The distribution of bar length peaks at values which are 
similar to those found by Marinova & Jogee (2007), Durbala 
et al. (2008) and Aguerri et al. (2009). However, the distri- 
bution presented by Barazza et al. (2008) peaks at a value 
which is about twice as low. This difference can likely be 
explained by differences in the samples and in the way bar 
length is defined. The sample used in Barazza et al. (2008) 
is at a redshift range of 0.01 < z < 0.03, while the sample 
studied here is at 0.02 < z < 0.07. Since they also used 
SDSS images, this means they have on average better spa- 
tial resolution to resolve shorter bars. In addition, they used 
the position of the peak in the ellipticity profile in the bar 
as a measure of bax length, which is known to usually result 
in lower values than other definitions of bar length (see e.g. 
discussion in Athanassoula & Misiriotis 2002; Erwin 2005). 
It is also interesting to note that no bar extends to a radius 
larger than « 3 times the disc scale length, and that almost 
all bars are shorter than r24 and R90. The distribution of 
the bar-to-total luminosity fraction shows a peak at « 0.1. 
The corresponding peak for the mass fraction is at a higher 
value, as expected, since bars are usually mostly made up of 
old stars. As one can see, bars can contain up to about 40% 
of the total galaxy stellar mass. 
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Figure 2. Distributions of bar properties for galaxies with classical and pseudo- bulges. Bars in galaxies hosting classical bulges tend 
to be larger than those in galaxies with pseudo-bulges. This is seen for all parameters used to define bar size, although the overlap is 
significant. The bar-to-total luminosity ratio is also larger in galaxies with classical bulges. Considering all other measured bar properties, 
including bar-to-total mass ratio, galaxies with classical and pseudo-bulges host similar bars. The arrows in the I/^ar panel indicate the 
median values of the corresponding distribution for classical (red arrow) and pseudo-bulges (blue arrow). 



3.2 Bars in galaxies with classical and 
pseudo-bulges 

As mentioned in the Introduction, classical bulges seem to 
be formed in relatively fast and violent processes, such as the 
merging of smaller units, as opposed to pseudo-bulges, which 
are thought to be formed from disc instabilities, such as bars. 
It is thus important to ask if bars in galaxies hosting pseudo- 
bulges are in any sense different from their counterparts in 
galaxies hosting classical bulges. This question is answered 
in Fig. 2. The only meaningful difference I find concerns bar 
length. Bars in galaxies hosting classical bulges tend to be 
longer than bars in galaxies with pseudo-bulges. This is con- 
sistent with previous results that show that bars are longer 
in early-type disc galaxies (e.g. Erwin 2005; Aguerri et al. 
2009), since such galaxies generally host classical bulges, al- 
though there are many examples of early-type disc galax- 
ies hosting pseudo-bulges (see e.g. Erwin et al. 2003; Lau- 
rikainen et al. 2007). There is also a trend in which bars in 
galaxies with classical bulges have larger bar-to-total lumi- 
nosity ratios, but this vanishes when one considers the mass 
ratio. 

One can also ask whether pseudo-bulges relate to their 



corresponding bars differently than classical bulges, and this 
is explored in Fig. 3. It shows that, for a fixed bar-to-total 
mass ratio, pseudo-bulges are on average clearly less mas- 
sive than classical bulges. Evidently, this relation arises from 
(i): as shown above, most bar properties are unrelated to 
bulge category, in particular Bar/T (mass), and (ii): the 
separation between bulge categories tend to pick less mas- 
sive bulges as pseudo-bulges, and, conversely, more massive 
bulges as classical bulges. Note, however, that pseudo-bulges 
selected this way also have significantly lower Sersic indices 
and higher star-forming activity, as compared to classical 
bulges (see Paper I). This indicates that the separation is 
natural, since such dichotomies are expected between classi- 
cal and pseudo-bulges. In other words, Fig. 3 is just another 
way of looking at these previous results. It stresses that if 
pseudo-bulges are formed through disc instabilities such as 
bars then indeed one needs more than disc instabilities alone 
to build a classical bulge, as bars in galaxies with classi- 
cal bulges do not seem more apt to originate more massive 
bulges than bars in galaxies with pseudo-bulges. 

Conversely, if bars are generally quite similar in galaxies 
with classical and pseudo-bulges, as Fig. 2 indeed indicates. 
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Figure 3. Stellar mass of the bar divided by stellar mass of the 
galaxy (i.e. the bar-to-total mass ratio) plotted against bulge 
mass, for classical and pseudo-bulges, as indicated. Clearly, for 
a fixed bar-to-total mass ratio, pseudo-bulges are on average less 
massive than classical bulges. 



then there is no reason to beheve that secular effects related 
to the evolution of the bar do not occur in galaxies with 
classical bulges. In fact, since bars in the latter tend to be 
longer, their effects might well be stronger. This suggests 
that the possibility of having composite bulges, i.e. galaxies 
hosting both a classical and a pseudo-bulge, is actually quite 
likely. This has been explored in Gadotti (2009b), where 
some bulges were found to have mixed properties: while they 
appear classical from a structural point of view, the presence 
of a spectral signature of young stellar populations is typical 
of pseudo-bulges. Such combination of features naturally ap- 
pears if one has a more massive and extended classical bulge 
component, in which a small and young pseudo-bulge is em- 
bedded. The classical bulge component produces an imprint 
in the structural analysis, while the young stellar population 
from the pseudo-bulge component shows up in the spectral 
analysis. It seems natural that a disc galaxy can undergo 
the processes that lead to the formation of a classical bulge 
and those that lead to a pseudo-bulge. Curiously, Nowak 
et al. (2010) describe two examples of galaxies with com- 
posite bulges, in which the classical component is actually 
smaller. 



3.3 Bar strength and growth 

In this section, I will look for correlations using the struc- 
tural properties of bars and their host galaxies that can 
serve as a test to the theoretical predictions on the secular 
growth of bars. Essentially, one should check if bars indeed 
grow longer and stronger in time, and thus one first needs 
to define bar strength. One readily accessible measure of 
bar strength is the bar ellipticity. All things equal, a more 
eccentric bar induces stronger non-axisymmetric forces and 
torques on the otherwise close to axisymmetric potential of 
a disc galaxy (see e.g. Martin 1995, and references therein). 
Thus, bar ellipticity is a direct measure of bar strength and 
does not depend on any other galaxy property. A more so- 
phisticated measure of bar strength involves directly esti- 
mating the torques induced (Block et al. 2004; Laurikainen 
et al. 2004; Buta et al. 2005). Such torques are normalised by 



the axisymmetric component of the galaxy potential, which 
means that, all things being equal, a bar in a galaxy with 
a more massive bulge is weaker (according to this estimate) 
than an identical bar in a galaxy with a less massive bulge. 
The same bar can thus have a different strength, according 
to this definition, depending on the galaxy it resides. The 
torque is usually measured at the point in which tangential 
forces along the bar are at a maximum. Alternatively, it can 
be measured at other specific points. For instance. Grouchy 
et al. (2010) relate rings and bars using bar forcing at the 
location of the rings, rather than the position of maximum 
bar torque. Thus, it can be said that this estimate may be 
more related to the impact of the bar on the overall galaxy 
evolution, rather than the strength of the bar itself, and it 
also depends on which point along the bar the torque is mea- 
sured, and that can be adapted depending on what is being 
investigated. Therefore, in the context of this study, it can 
be argued that bar ellipticity is a more useful, global mea- 
sure. This is especially true as the theoretical prediction to 
be tested is the evolution of bar strength along the history 
of individual galaxies, and the axisymmetric component in 
the potential of a given (isolated) galaxy is not likely to 
change in time as much as the difference in such component 
for different galaxies. 

Interestingly, however, bar torques are very well cor- 
related with bar ellipticity (Laurikainen et al. 2002; Block 
et al. 2004), and therefore the results presented below are 
not likely to depend significantly on how bar strength is eval- 
uated. It should also be noted that a recent study (Comeron 
et al. 2009) suggests that a combination of both bar torque 
and ellipticity can be a useful measure of bar strength. 
These authors found that the curvature of dust lanes along 
bars, which is also related to bar strength (see Athanassoula 
1992) , is better related with such combination than with bar 
torque alone. 

Evidently, bar mass, length and even bar boxiness are 
related to bar strength. In fact, the simulations in Athanas- 
soula & Misiriotis (2002) indicate that bars grow stronger in 
time by getting longer, more eccentric and more boxy, which 
is consistent with the results from Fig. 4, which shows that 
bar ellipticity and boxiness are correlated (for both galax- 
ies with classical and pseudo-bulges). Since both parameters 
are related to bar strength, I will use the product of the two, 
i.e., e X c, as a proxy. It is important to note that the trends 
discussed below are stronger when one uses ex c rather than 
either e or c alone. However, there is no improvement if one 
uses e X c times the bar mass or the bar-to-total luminos- 
ity and mass ratios. Thus, the choice of using the product 
e X c as a relevant parameter related to bar strength is not 
arbitrary, but based on the fact the results become clearer 
by taking the boxiness c into account, as opposed to other 
parameters. A comparison between e x c and bar torque is 
beyond the scope of the present paper and will be published 
elsewhere. 

I now need a working definition for bar size. Figure 5 
shows, consistently with previous results, that bar size cor- 
relates with galaxy size, and this is regardless of whether one 
uses the bar semi-major axis Ltar, or the bar effective radius 
re,bar, for the former and h, r2i or R90 for the latter. It also 
shows that Ltar and r^^i^ar are well correlated, as expected. 
Because larger bars are in larger galaxies, a definition for 
bar size in the context of this study has to be normalised by 
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^ar (^P'^) R90 (kpc) 

Figure 5. Correlations between the diflferent parameters used to define bar size. They show that, as expected, larger bars reside in larger 
galaxies. In addition, bar eflfective radius is equivalent to bar length, and h, r2A and i?90 can be similarly used to obtain normalised bar 
sizes. 
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Figure 4. Correlation between bar ellipticity and boxiness for 
galaxies with classical and pseudo-bulges, as indicated. 



galaxy size. Figure 5 shows that one can use any combina- 
tion between Lbar or re,bar and h, r2A or R90 to define the 
normalised bar size. 

Figure 6 shows that there is a positive trend between 
all applicable definitions of normalised bar size and bar e x 



c for both galaxies with classical and pseudo-bulges. This 
is consistent with the expectations from theoretical work, 
although clearly there is considerable scatter in some of the 
plots. 

Figure 7 shows that bar size also correlates with bulge 
size (see also Athanassoula & Martinet 1980), which is per- 
haps not too surprising, but, interestingly, this correlation 
depends on bar ellipticity. In fact, the correlation is weaker 
for bars with e > 0.7. Surprisingly, re/re,bar shifts to lower 
values for more eccentric bars. This effect is consistent with 
the previous finding that more eccentric, stronger bars are 
longer, and will be further explored below. 

We have already seen that longer bars have higher ex c 
values, but the theoretical work I aim to test here also pre- 
dicts that such bars are dynamically old, i.e. they had more 
time to evolve. This would happen, for instance, to a bar 
that resides in a galaxy that have reached a dynamically ma- 
tured state earlier.^ Given the current cosmological picture 
of galaxy evolution, more massive galaxies reach this point 
earlier (something referred to as "downsizing" - see discus- 
sion in Sheth et al. 2008, see also Cameron et al. 2010). But 
more massive galaxies tend to be those in which the bulge- 
to-total ratio is larger (see e.g. Gadotti 2009b), and thus if 
the theoretical work is correct - and if one assumes that 
B/T indicates bar dynamical age - one should find a trend 



* This assumes that bars are long-lived structures and not recur- 
rent as in the picture proposed by Bournaud & Combes (2002). 
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Figure 6. All dlflferent measures of normalised bar size plotted against bar e x c for galaxies with classical and pseudo-bulges, as indicated. 
Although a few plots display large scatter, the positive trend between the parameters is clear. 
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Figure 7. Histograms of the ratio between bulge effective radius 
and bar effective radius for bars in three bins of ellipticity, as 
indicated. It is clear that both parameters are correlated, at least 
when e < 0.7. Note that the peaJjs in these distributions shift 
successively to lower values from bars with e < 0.6 to bars with 
0.6 < e < 0.7 and finally to bars with e > 0.7. 



in the sense that longer bars are found in galaxies with larger 
B/T. This is consistent with previous observational results 
that longer bars tend to be in early- type galaxies, but since 
I have measurements of B/T I can test it directly. 

Figure 8 shows the different measures of normalised bar 
size plotted against bulge-to-total ratio. This is done for bars 
in three bins of ellipticity, as in Fig. 7, for reasons that will 
be clear shortly below. Although there is considerable scat- 
ter in some of these plots, one sees clearly that longer bars 



tend to be hosted by galaxies with more conspicuous bulges, 
which is thus consistent with the theoretical expectations. 
If these plots are done with all bars grouped together, re- 
gardless of their cUipticities, this relation is more difficult to 
see. The reason for that is that bars with different elliptici- 
ties follow somewhat parallel lines in these plots, and their 
combined scatter dilutes the relation. This is easier to see 
in the bottom right panels in Fig. 8. The top panel shows 
LBar/r24 plotted against B/T, with colour coding indicat- 
ing bar ellipticity. The black solid line is a fit to all points. 
Most green points (bars with e > 0.7) are above this line, 
whereas most red points (bars with e < 0.6) are below this 
line. Black points, representing bars with 0.6 < e < 0.7, lie 
along the line. The bottom panel shows separate fits to these 
data points, separated in different ellipticity bins. One sees 
that the fits describe roughly parallel relations. It should be 
noted that this effect is apparent for all measures of nor- 
malised bar size. 



4 DISCUSSION 
4.1 Possible bieises 

Some of the physical parameters explored in this study are 
difficult to measure, such as bar boxincss, in particular con- 
sidering that the physical spatial resolution of the images 
used is relatively poor. In Gadotti (2008), I presented a 
study using images of 14 nearby gala:xies {z 0.005 - 12 
of them barred galaxies) and images of the same galaxies 
artificially redshifted to z = 0.05, which results in the same 
physical spatial resolution of the sample used here. The red- 
shifted and original images were decomposed using the same 
methodology as in Paper I. By comparing the results from 
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Figure 8. All different measures of normalised bar size plotted against bulge-to-total ratio, for bars in three bins of ellipticity, as 
indicated. Although some panels show considerable spread, in most of them a clear relation is seen in the sense that longer bars tend 
to reside in galaxies with more conspicuous bulges. Bars with different ellipticities describe roughly parallel lines in this relation. This 
is better seen in the lower right panels: the top one shows the data for all bars in the panels where is plotted against B/T, 

colour coded by bar ellipticity (the solid line is a fit to all points); the bottom one shows separate fits in the different ellipticity bins (the 
dotted line in the bottom panel is also a fit to all points, plotted again for reference). 



such decompositions, it was verified that parameters such as 
disc scale length, bulge effective radius, bulge-to-total ratio 
and bar-to-total ratio can be reliably estimated in the low 
resolution regime. Bulge parameters are particularly sensi- 
tive to seeing effects, though, if < 0.8 x iPWHM, but, as 



shown in Paper I, 97% of the bulges in the sample have re 
above this threshold. 

I can now use the results from the decompositions in 
Gadotti (2008) to verify if the other bar structural parame- 
ters discussed here can also be reliably obtained in this low 
resolution regime. Figure 9 thus compares the estimates for 
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bar Scrsic index, length, ellipticity and boxiness from the 
original and artificially redshifted images of the 12 barred 
galaxies in Gadotti (2008).* It shows that estimates of such 
parameters are indeed robust even at the resolution of the 
images used here. No clear systematic biases arc seen, ex- 
cept for a tendency for lower ellipticities, as measured in the 
redshifted images. This likely results from the rounding and 
dilution of the bar light at low resolution. It might suggest 
that the peak in the distribution of bar ellipticity in Fig. 1 
is slightly offset to lower values. However, the difference is 
on average only a few percent, and is the same through the 
whole range of values obtained, and thus not harmful to the 
results presented here. 

One could also be worried that bars might show up more 
eccentric and boxy in galaxies with less conspicuous bulges, 
where there is less contamination from bulge light in the bar 
isophotes. While this is a concern for results based on ellipse 
fits, it should be stressed that the bulge light is modeled in 
image decomposition, and thus such concern does not apply 
here (see fig. 6 and discussion in Gadotti 2008). Further, the 
results shown above point out the opposite: more eccentric 
and boxy bars are longer and hosted by bulge-dominated 
galaxies. 

Finally, in the same context, one could argue that higher 
values of ellipticity and boxiness are easier to obtain in the 
case of longer bars, since these have a larger fraction of 
themselves away from the bulge dominant light. However, 
Fig. 9 shows that boxiness can be reliably measured at the 
low resolution regime, i.e. even when bar light is diluted. 
Furthermore, considering the four galaxies that have their 
ellipticities underestimated at the low resolution regime (see 
lower left panel in Fig. 9), one sees no tendency for them to 
be particularly short (these galaxies are NGC 4314, 4394, 
4477 and 5701). 



4.2 Bar length in simulations 

As discussed in Erwin (2005), to compare bar lengths in ob- 
servations and simulations one should use LBa-r/h. He pro- 
vides this value for a number of models from the literature 
(his Table 8) and argues that there is clearly a lack of short 
bars in simulations, as compared to observations. Consid- 
ering only those models in which a live dark matter halo 
is used (Berentzen et al. 1998; Athanassoula & Misiriotis 
2002; Valenzuela & Klypin 2003; HoUey-Bockelmann et al. 
2005) one sees that indeed most of them produce bars which 
are in high end tail of the distribution of I/Bar/ft shown in 
Fig. 1. In fact, Berentzen et al. (1998) and Athanassoula & 
Misiriotis (2002) report results in which Lear/Zi > 3, a value 
not found in the sample studied here, which is nevertheless 
biased towards large bars. 

The models that produce the shortest bars are models 
Al (LBar//i = 1.4 - 1.5), A2 (LBar//i = 1.1 - 1.3) and B 
(LBar//i = 0.8 - 1.0) in Valenzuela & Klypin (2003). The 
main difference between models Ai and A2 is that the latter 
is kinematically colder than the former in the central region. 
Model B has a smaller halo and a more massive disc. It 



* Note that two bars in this sample become overly faint in the 
redshifted images and axe thus not fitted in these images. 
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Figure 9. Bar structural parameters, as determined witli the 
redshifted images, plotted against the same jjarameters obtained 
with the original images. The dashed lines indicate a perfect cor- 
respondence. I/Bar is in axcseconds, and measurements from the 
redshifted images are scaled back to the original galaxy distance. 
No clear systematic biases are seen, except for a trend in which bar 
ellipticities, as measured in the redshifted images, have slightly 
lower values, compared to those measured with the original im- 
ages. The difference is on average only a few percent, and is the 
same through the whole range of values obtained, and thus is not 
harmful to the results presented here. 



appears that a less efficient transfer of angular momentum 
between disc and halo is the cause of these bars being short. 

A related issue is that at least in some models (see e.g. 
Martinez- Valpuesta et al. 2006) long bars are formed quickly, 
i.e. after 1 — 2 Gyr only, while the results above suggest that 
long bars arc old. As both bar size and disc size change in 
time, it would be most useful if theoretical work presents 
how L^ai/h (or LBar/T'24) evolve with time. 



4.3 Secular evolution of haxs 

Although there seems to be some quantitative disagreement 
between theory and observations concerning bar sizes, a 
plausible interpretation of the results above suggest a qual- 
itative agreement. The data show that longer bars (using 
normalised size) have higher values of e x c (Fig. 6), and 
reside in galaxies with higher B /T (Fig. 8). Now, one can 
interpret these results, assuming that higher values of B/T 
indicate dynamically older bars. This would thus suggest 
a general trend showing bars getting longer and stronger 
with time, generally consistent with sirrmlations. This is not 
the whole picture, though. As discussed in the Introduction, 
some simulations show that bars can get very large quite 
fast, quickly become shorter at the onset of the first verti- 
cal buckling instability, and then start growing slowly again. 
Furthermore, simulations taking the effects of gas into ac- 
count show that the secular growth of bars can be halted, 
or even reversed, if the gas content is enough. 

Since the sample studied here is biased towards large 
bars, it is plausible that many of these bars have evolved 
enough, and already went through the first vertical buckling 
instability. For the same reason, and also because all galax- 
ies in the sample have stellar masses above 10^° M©, the gas 
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content in these barred galaxies is likely on the lower end of 
the corresponding distribution. This might have helped to 
make the trends in Figs. 6, 7 and 8 more apparent. Never- 
theless, the complicating effects from gas dynamics and the 
vertical buckling are probably contributing to the spread 
seen in these figures. A similar work, with a sample of short 
bars, in gas- rich system, is likely to shed light on these issues. 

With a different methodology, Elmegreen et al. (2007) 
found that normalised bar size correlates with bar strength 
and galaxy central density, and also conclude that their re- 
sults "suggest that bars grow in length and amplitude over 
a Hubble time" . 

Observational results like those in e.g. Rautiaincn et al. 
(2005, and references therein) indicate that most bars end 
near their corotation radius. Along with our current under- 
standing of the orbital structure in baxred galaxies, this tells 
us that bars cannot grow longer if they do not slow down. 
Therefore, if the interpretation of the results discussed here 
is correct, namely that bars grow longer with time, one can 
speculate one step further and suggest that bars also slow 
down with time (which is again qualitatively consistent with 
theory), and with bulge prominence. Such a relation between 
the bar pattern speed fie and B/T would be in the sense that 
bars rotate slower in galaxies with more prominent bulges, 
since these galaxies have longer bars. This relation could 
translate to a dependence of fls on Hubble type, although 
there is some scatter in the relation between B/T and Hub- 
ble type (e.g. Laurikainen et al. 2007; Graham & Worley 
2008). On direct observational grounds, no solid conclusion 
can currently be drawn about a dependence of bar pattern 
speed on Hubble type, and this is due to the difficulty of ob- 
taining reliable estimates of fis, in particular for late-type 
galaxies (see e.g. Gerssen et al. 2003; Treuthardt et al. 2007; 
Corsini 2010). 

An outstanding and unforeseen new result in this study 
is the existence of roughly parallel tracks in the correlation 
between bar normalised size and B/T, for bars with different 
ellipticitics (see Fig. 8). A straightforward way of interpret- 
ing the existence of these parallel tracks is to conceive that 
bars could form with different normalised sizes and elliptic- 
itics, and then follow a somewhat parallel growth. This is a 
new aspect that can be investigated with theoretical work. 
It would be very interesting to study bars formed with dif- 
ferent normalised sizes and ellipticities in simulations, and 
check whether they follow similar evolutionary paths. 



5 CONCLUSIONS 

I have explored the results from detailed 2D image decom- 
position of nearly 300 barred galaxies with stellar masses 
above 10^° Mq, at « - 0, concerning mainly the structural 
properties of bars. This results in a thorough description of 
bonafide stellar bars in the local universe, including distri- 
butions of bar Sersic index, ellipticity, boxincss, length and 
bar-to-total luminosity and mass ratios. The interplay be- 
tween bars and the bulges and discs in their host galaxies 
was also examined. Such detailed characterization of local 
bars can be compared with similar results from studies with 
samples at higher redshifts, in order to directly investigate 
how the structure of barred galaxies evolve in time. It can 
also be compared to results from theoretical work, in order 



to assess how well theory can describe the structural prop- 
erties of such bars. Furthermore, a complete description of 
the properties of local bars is useful to insert ad hoc models 
of bars in a theoretical framework. 

Bars in galaxies hosting classical and pseudo-bulges 
share similar properties, except that bars in the former are 
on average larger than those in the latter, considering both 
absolute and normalised sizes. This is consistent with previ- 
ous results comparing bar sizes in early- and late-type disc 
galaxies, as usually the former host classical bulges, whereas 
the latter host pseudo- bulges. For a fixed bar-to-total mass 
ratio, pseudo-bulges are on average significantly less massive 
than classical bulges. This indicates that, if pseudo- bulges 
are formed through disc instabilities such as bars, then more 
than that is necessary to build a classical bulge. 

Normalised bar size is correlated with the product of 
the bar ellipticity and boxiness (e x c, which is related to 
bar strength) and with bulge-to-total ratio, B/T. If higher 
B/T indicates more evolved bars, then these results can be 
interpreted as qualitatively consistent with general expec- 
tations from theoretical work, which suggests that evolved 
bars grow longer and stronger in time. This would come 
along with a decrease in bar pattern speed, and assumes 
that effects caused by gas dynamics are sufficiently small. 
Bars with different ellipticities follow parallel tracks in the 
trend found between normalised bar size and B/T, suggest- 
ing that bars could form with different normalised sizes and 
ellipticities but still follow similar evolutionary paths. 
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